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ABSTRACT 


Galea O72 ye Arctic numerical model was used to 
oo Oremunc  stonmiiticant dynamics of the Bering Sea. The 
mode | wasesetsupeand tested for numerical stability. 
Climatological wind fields were used along with source- 
Sink exchange transports to drive the model. The runs 
investigated the results of adding or deleting bottom 
Lieionmenon-linecarities and bathymetry. The results 
Heine aepy etic model were compared to a circulation 
tieremmmpueposcd by Hughes (1972). The model showed 
that the circulation is strongly bathymetry-dependent and 


primarily driven by the sources and sinks. 
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I. INTRODUCTION 


ie mOMycet@eOLtmElIs Proj;ject was to conduct a numerical 
expueoration Of large-scale flow ofethe deep Bering Sea 
with an initial goal of developing a working numerical 
ieee ror the desired area, Following this, the most 
important goal was to study the effects brought about by 
varying the important oceanographic ET Dili s an themmodel. 

iMemea@sresmodel utilized was one developed by Galt 
ie ees 7>) fOr use in the Arctic Ocean. This model 
assumes homogeneous water, has capabilities of varying the 
depth, and includes lateral and bottom friction. The flow 
in the model is driven by a wind stress applied at the 
Surface and the introduction of sources and sinks into the 
vorticity field at the boundaries. This model then was 
ani ounmete1or Uwse in the Bering Sea by alteration of the 
grid field to conform with the study area. 

Figure 1 shows the entire Bering Sea basin and shelf 
area. Two main factors influenced the reasoning used in 
choosing the deep water basin part of the Bering Sea area 
for the model. The first was due to a combination of the 
bathymetry of the Bering Sea and the expected physical 
Himitations of the model. Beyond the northeastern boundary 
Cemenecmnode! the Bering Sea bottom 1S a fairly long, flat 
continental shelf which presumably will not drastically 


aitect the dynamics of the western portion of the Bering 
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coca. Due to the physical configuration of the model area, 
the mathematical constraints and geometrical limitations 
placed upon the model [Willems, 1972; Galt, 1973], a larger 
area was avoided. The second reason deals more with the 
Peele Ol senewoasSilemprO}ECL, Existing obserwed and 
MapoLicmmenml data do not extend much further beyond the 
edge of the continental shelf than the model does and, 
considering that the model flow was to be Peinieed to 

this data, the model study area was determined as shown 
(ieemre 2). 

Mice Ince seri1d System wsed 1n the model was based on 
iMaomanmeutar cell as used by Galt (1972). The resolution 
eretiemerid was 1/2 of latitude at 58° North latitude or 
Pporoxamately 30 nautical miles (55.6 kilometers). In 
the model, the grid was coded as a parallelogram with N 
columns and M rows, thus giving a total of N x M points 
(Z4 x 285— 672). the grid points were numbered sequen- 
Piolly tromelett to right by rows. Two additional vector 
Pens Weremurilazed tO Specify the extent of the interior 
domain aleamene dimensioned M. The ees aemmaiiy studied 
was coded by use of initial boundary conditions inside the 
Original parallelogram. 

AS PiewtoOusly stated, one of the primary driving forces 
for Galt's model is the application of wind on the surface 
aero nl icmremial Wind data utilized for this inves- 


Mit MovmeNacmobtained=:rom Aagaard [1975]. Source-sink 
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betma@aryeeonditions were taken from Hughes [1972], 
hereafter referred to as Hughes, and Coachman [1973]. 

The Bering Sea has been studied for years but, until 
mie work Of Hughes, a detailed study of the entire water 
~Epollatitom pattern based on observed data had not been 
mone. Onemwos the results of the several studies conducted 
wii wind anppears to be the primary driving force 
mea the sliriace water circulation in the Bering Sieg. iste: 
CimviiiCmerclas UuSed for this project were obtained as 
monthly climatological means taken from atmospheric pres- 
Sea omiiavce ss lhe obvious point of interest in these fields 
Veomene strong increase in force in the winter months. The 
basic directional driving force, however, remained the same. 
This maintained a generally cyclonic flow pattern present 
Mime icmourrace circulation, 

There had been three basic methods for the determination 
Sumetrculation im the Bering Sea prior to Hughes' work. 

The first was the classical dynamical method used for 
GCOmputing geostrophic circulation. This method has never 
moo llivwerpmoauced jonelusivye results due te tneescattered 
and sparse data field available at any one time. The 
second method has been by water mass analysis. Neither 
eore nor volumetric methods have been used and the isen- 
tropic analysis has received little use. The waters of 
Bie Morim@nocasare Gharacterized by a cold intermediate 
Pyerepumenes| ana @ samilarity to the North Pacific 


Wem Sermreminetwe deeper layers [Sverdrup, 1972]. This 
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then leads to the difficulty in analysis. The third 
method has been the use of wind-driven models. These 
have been difficult to establish as a forecasting tool 
due to the rather inadequate data available and the 
difficulty in the handling of the averaging methods used 
on the atmospheric pressure maps [Hughes]. 

Vanouemeanpamese, U.S.S.R., and U.S. investigators 
have looked at separate sections of the aie pattern in 
the Bering Sea over the years and, as a result, extensive 
data files have been gathered by these countries. Hughes 
combined parts of this data with three years of his own 
direct observations and measurements, and conducted a 
Grercuem study of the single-layer circulation of the 
Bering Sea deep-water basis. As his measurements were 
all made during the mid-summer months, Hughes directed 
tseattention mainly to the flow during that period. 
Based on other data and theoretical work, Hughes does 
Speculate on the flow during the winter months but draws 
[mows pecdmete Conclusions. 

Hughes! work was principally the result of drogue 
studies and in-situ measurements. Secondarily, Hughes 
compiled transport flow values through the various straits 
and around the boundary areas in sufficient detail that a 
direct application was made to the model where necessary 
for sources and sinks. To make these calculations, Hughes 


had to make a rather exact water budget. He obtained 
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cross-sectional areas of the straits and passes which were 
considered major inflow and outflow areas. He then calcu- 
Haeca Water transports and compiled the final budget 

emi ezinee tie concept of continuity. 

Hughes coupled his primary study with a qualitative 
emelysis Ot bathymetric effects on the circulation. He 
Showed that, due to the Shirshov and Bowers Ridges (Figure 1), 
the basic cyclonic flow was, in fact, broken -:into several 
ditferent gyre patterns surrounding the ridges. Even so, 
he found that the general overall cyclonic pattern remained. 

With Hughes! work as a starting point in the under- 
standing of the flow system in the Bering Sea, the model 
Was then used to explore the significant dynamics in the 
deep water basin. The results of the model] runs were com- 
pared with Hughes' work on a qualitative basis only. 

There were three fields of variables that required 
investigation before the model could be run efficiently. 
Firstly, there were numerical constraints mach Peer te 
variable items that control the basic model program and 
Smevireietinal flow pattern results do not depend. Such 
an item would be the time step used. Secondly, there 
Were cCONLipurational constraints, which in this case 
Peretrecde tothe geometry of the model. Such items as grid 
Tiiem sates in this group. The third group consisted of 
the program parameters. As examples, this list included 
hicevamrabte depths at cach grid point, the variable 


fidecional terms, and the variable driving functions. 


ie 





The variable parameter depth is coded as a depth factor 
WimemecanmuesVaried to include the entire topographic pro- 
MeomOormcaiy percentage thereor, iancluding a flat bottom case 
SPeezeromperceim. | [his factor alone offered an opportunity 
fomcexplore Hughes’ analysis of bathymetric effects. At the 
Same time as this was being accomplished, the effects of 
Vamameetatenal and bottom friction and the inclusion and 
PMemusi@lvOL 1On-linearitics was explored. Exploratory 
runs were also conducted to determine the most efficient 
values of the numerical and configurational constraints to 
be used. 

Sensitivity tests were made on the model to determine 
DeseStability with the different parameter field variables 
working in the model. Throughout these runs, a hypothetical 
wind force was used as the driving force. Following these 
Picweamcno1ce was made of the best runs and these were put 
together with a forcing function derived from climato- 
Hocvedmewima data to form the final series of mune: These 
then were compared: with Hughes best prognosis of the overall 
circulation. The next section will deal with the equations 


used and the actual set up of the model. 
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a EMAL ICS 


The area cnosen for study includes all of the major 
topopraphical features encountered in the Bering Sea; 
foncanental snelf, ridges, deep basins, straits, and closed 
and open boundaries. As has been stated, the major driving 
forces in the model were wind stress on the surface and a 
source-sink exchange along the boundaries. Hughes concluded 
that although the relative magnitudes of the wind stress 
components were much less in summer than in winter and that 
even in some cases the wind field became anticyclonic, the 
Sveralim circulation ot the Bering Sea maintained a cyclonic 
Daeecrnaws ice also Showed that the source-sink transport 
exchange between the Bering Sea and the Pacific Ocean for 
the summer months was significant primarily only in the 
Kamchatka and Near Straits and somewhat in the Buldir Pass. 
Geewas sugeested that exchanges along the other boundaries 
were small and of Secondary significance in the overall flow 
scheme for the summer months. 

inemcucepewaver basin Of the Bering=scea is divided by 
the Shirshov and Bowers Ridges (Figure 1). The depth of 
water above the ridges is less than 1000 meters and in the 
basins to either side of the ridges is greater than 3000 
meters. Depth within the model was coded as a deviation 
from the average depth (taken in this case as 3000 meters). 


eee moi vaimess were Normalized to this average depth 
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before entering them in the program. An algorithm was 
then used that allowed the model to work with percentages 
Srernewnonmalized deviation factors. For example, speci- 
mallee zenOepercent resulted in all points having a depth 
value of one and the model treated the basin as a flat 
bottom case. 

The model used was originally developed as a homo- 
geneous, barotropic model because a baroclinic model 
would involve much greater complexities in the mathematics 
than was desired at this stage. In addition, Galt (1972) 
has shown that there may be depth independent components 
to the current, and thus, this seemed a useful starting 
PeiiestoOr a numerical investigation. 

The frictional terms of the equations used in the model 
PaowOLmuvoOmuasicG types, One 15 lateral friction and the 
Senereisebpottom friction. Although both terms appear in 
prcecduatmen, they are both truly effective only in parti- 
lar anedame sine lateral friction term effects are noticeable 
along the boundaries when no slip conditions are specified. 
imma tree Slip condition there is no Doundary for the lateral 
PoroerTone te interact with. Bottom friction is dependent upon 
and inversely proportional to the depth and basically deals 
TEI MEnemcmSSipation Of VOrticity. Because of its relation 
Pomc ieomuepi i. it lS Cttective primarily in the relatively 
Shallow areas. It should be noted here that the minimum 
feet teizedmin the model 15 0.01 (30 meters) of the 


normalized averaged depth. 
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The following development is intended to give the 
Wem@enma CUNSOry inspection of the basic equations and the 
methods used to manipulate them. For a more detailed 
Mitioti oe aeroumeme reader 1S referred to Galt (1972, 1973). 
The basic integrated form of the equation of motion was 


mre Starting point, 


eens, hvu AN {1} 
Dv 4 gue —LOP 4 Kgtv— Bu 4 Te (2) 
Ot Pos h ph 

where: 


u and v are horizontal components of velocity indend- 

Epiemon depen; 

Pis density and is constant; 

+* and T? are wind stress components; 

leeseitane depth; 

f is the coriolis Panameter,; 

soem and £ are functions of position; 

K and R are constants specifying the effectiveness of 

hoOrdaoneale and svertical trictional forces respectively. 
Along with these two basic equations of motion, the continuity 


equation, 


Sx bul + $5 lhy}=0 (3) 


Peele transpOmt Stream function equations, 


_ av _ av 
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Vetescembinea With some manipulation to form the vorticity 


equation, 


af _(9 VR} 0 (EEE) note 8 e+ry-o())| 





where: 
a 
g= Wu 
ol) —9 
Vi Jeigg ty a 


imei ceoeinitial and boundary conditions established 
in the beginning of the program, the preceding equations were 
Pouca erorevOrticity and Stream function. The three conditions 
which must be specified are: 

a) WU -given within the region of interest at t=0; 

b) W -given oa the boundary of the region for all 

; time, (this one is the same as establishing 

Seunce-Silk bOlndary conditions) ; 

G) € -given Wit hineene region Ofsanterest at t=0. 

Reiter wem-dimensionalizing all quantities in the fore- 
going expressions, an application of the Adams-Bashforth 
PimeeoMevel! mminite differencing scheme was used to integrate 
the equations. The following non-dimensionalized form of (5) 
was used: a§ Bee), 

a = 


1§ 





a 


which was integrated with the following scheme: 


g(r rar l=St}+ ea t-sa(f-aryjor (7) 
In these last two expressions g' is the finite difference 
monn Of : 
=(yxw'}-v(28t4| + pv'e-Zlervy' ute i 
+7 «(F| (8) 


The expansion of the terms in (7) contain three non-dimensional 
SoetrycTenesewnich povern the scaling of the solution and 


dictate the relative importance of various terms. They are: 


W 


l= SDF ma lenmeovenms Mon-linearitics. Jf Cl = 0, 


the model considers a linear problem only; Q is the 
— el 
= =, Cea Pavey olen ie t 1 On; 1s the 
= (3 Buc e 6) 


oe. Ekman number; 


=y=t. EOMmmcols DOLtTOM Pricer ron, my is the bottom 
drag coefficient. 

Where ; | 
[ei omEnc Mormalyzed COrlOl1S parameter utilizing the 
North pole as a base; 
nero tne finite Pecumcnce men (55.0. Kil. ) ; 
D is the average depth of the model; 
We Pmecdwualeroul overdrup, (0° Wo /eoe. 
iieecmisealiso a fourth term in equation (7) which 1s of 


importance. This term represents the torque added per unit 
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time by the wind stress. This 1s assumed constant in time 
and calculated only once at the beginning of the program. 

The original integrated Sverdrup transport data used in 
the actual runs was produced at the University of Washington, 
marcomptlingsclimatological pressure data for a period of 
Six years. This was converted to wind data and finally to 
Wind stress, all done manually (Aagaard 1973). The final 
form of the data, as it was used, was obtained by integrating 
in the X-direction from east to west as demonstrated by 
Sverdrup (1963). 

rat 1S: 


where: 

Hee Pema SS tran Soort inthe Y~direction; 

oeacetiie Wind stress; 

Bb is the beta function -2rcosg sl, 

9 iomeee wlatia tide ; 

PeTomehe radius of the carth; 

eicmecamel to 277 radians per day. 
Differentiating results in an expression for the curl of 


the wind stress, 1.e., 
curlT=My:f (1 O} 


This was the desired input for the model. An algorithm was 
developed which gave a finite-difference value of the inte- 


Peete oseranupetransport program inputs. Multiplying this 
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YATE)’ where H(t] Prcveomaerin dt miememmen point, curl 


of T was obtained. The only problem then remaining was 
the non-dimensionalizing of these terms in order to remain 
PoMolosteChemwiem the remainder of the program. This was 


aoe emul taplication factor of ——. 


F 


After integrating the vorticity equation foreward one 
muee step anl Cxpression relating vorticity and stream 
mone Lion similar to equation (6) was solved by a successive 
Gverrelaxation technique as described by Galt (1972). 
The optimum value of the relaxation variable used in the 
Neoeelewas saetcimined trom test runs. This configurational 
foro ines dependent on the grid size and number of 
grid points, satisfying the expectation that the relaxation 
Paaerunethon or the geometry of the model. 
Veronmeoneexecuting the series of runs to describe 
Me ebicotropr1c flow driven by the realistic wind data, 
it was necessary to run a series of exploratory tests on the 
Pimemicalwauc COntigurational constraints and the parameter 
field. The following section will deal with these exploratory 


runs and the sensitivity and stability tests made on the model. 
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III. EXPLORATION 





Once the model was put together in a running mode, 
a series of test runs were made in order to determine 
appropriate values for model variables and parameters. 
The first investigated was the relaxation constraint. 
Initially runs were made varying the relaxation value 
imoietec0eco 2.00, A plot of the relaxation constraint 
versus number of iterations per time step was made. A 
value of 1.80 was determined to maximize the convergence 
Gace tor the given model geometry. 

Wowie scmparically determined from Arctic model 
tests (Galt 1973) was initially used for C2, the con- 
Miomuincmoanametcr tor lateral friction terms in the 
equations just discussed. This value was derived uti- 
lizing a K value of 108 (in non-dimensional units), which 
is consistent with the lateral friction used in other model 
Studies. Throughout the test runs and later during the 
actual runs the value was varied periodically to check the 
model's sensitivity to this parameter. The different values 
used were on the order of one magnitude either side of the 
original derived value for C2. In all cases the model 
results, using these different values of C2, were interpreted 
as having less realistic energy buildup. Thus, the original 
value determined was used for all further runs. 


The time step value was investigated next. An original 
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value of 0.5 was chosen as a normalized value corresponding 
to a half pendulum day. Since F was scaled on the polar 
value (latitude = 90°), this was equivalent to taking a 
time step of twelve hours. Runs were made varying this 
numerical constraint to see what dependence, if any, there 
Was on the numerical stability of the model. When running 
ferme tne cul] bathymetry factor in the nee non-linearities 
Meanbotutom friction terms set to zero, and using a constant 
torque applied as a driving force, the model achieved 
Steady state solutions only with a time step of 0.0125 or 
less. With time-step values above this, including 0.5, 
numerical instability resulted. However, once the bathy- 
metry factor was lowered from 1.0 or 100%, the original 
emfemor 025 fave Satisfactory results in all cases tested. 
The value of 0.5 was chosen then as the standard time step 
to use throughout the model. 

DCMU Imnepilarities noted in stability at a 
depth factor of 1.0 in the last test series, the next 
examination was conducted to determine the model dependency 
werne batiymetry. In Galt*s (1972) original work, the flow 
in the Arctic model was relatively sensitive to bathymetry, 
particularly in the area of ridges. Hughes also showed 
mehitatuyely that the flow of the Bering Sea is affected 
by topographic relief on the ocean bottom. The kinetic 
energy, aS an output value in the program, was determined 
by summing the Pinertiecnwemcrey calculatedgat each grid point. 


i Was Considered to be a good value to monitor for detecting 
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Stabtiiteveor Mstability. In all cases, total kinetic 

energy was therefore monitored and proved to be a sensitive 
indicator of stability and displayed a dependence on bathy- 
metric modifications. Based on Hughes! work, it was expected 
that there would be a bathymetric dependency on the flow 
Ofethe Bering Sea. It was expected that steady state con- 
ditions would be achieved under approximately the same 
fomele1Ons aS in the Arctic. 

(ies irsteplot of kinetic energy versus bathymetry 
percentage was made on the model with Cl (non-linearities) 
omameoe(DGttonm £Y1iction) equal to zero and a constant 
torque driving the flow. In all runs where constant 
torque was applied, a force per unit area approximately 
Beewoyiescewas applied across the entire basin area. No 
Semuecacilk GOnaltions were applied in these cases. The 
Same runs were then redone with source-sink conditions 
applied along the western boundary of the study area, 
Specifically in the Kamchatka, Near and Buldir Pass Straits. 
From Hughes’ work, -the source-sink values applied were 
epeeltiedvas follows (Figure 3): 25 x 10° m”/sec. outflow 
through the Kamchatka Straits; 5 x i0° m/sec. inflow 
between the Komandorski Islands; 5 x 10° m/sec. outflow 
ipeememwestern section of the Near Straits; 30 x 10° eee. 
iMiehOowesSpheaa over the eastern portion of the Near Straits; 
echo X iol’ m/sec. outflow through the Buldir Pass. Figure 4 


ievomnoCoultsceot these runs. lt should be noted that the 


Detective niit to the data iS approximately exponential. 
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A determination was made on the basis of these plots that 
runs utilizing 20% bathymetry would satisfactorily describe 
the relevant dynamics. As a check on this decision, some 
footer tem=mugms were made at 505 bathymetry later on and 
no appreciable qualitative difference in the flow was noted. 

Plots were then made of kinetic energy versus normalized 
Piinicwoteps or the Same Set of runs. Initially all runs were 
terminated after 600 time steps, regardless of flow con- 
ditions. At that stage, the program produces a data file 
leielowotomed LOr tuture use,” This data file contains 
MiemextceIng Tun parameters at time step 600 so that runs 
Mays Se Started at that point and continued further Aa ‘Giner 

itesceplLotscewerenuncn extended to 1800 time steps (Figure 5). 
A stability determination was made from this type of 
Pot wrUmet@ac showed energy bulld-ups leveling off 

and reaching or, at least, definitely approaching asymptoti- 
cally a steady state condition with increasing time was 
considered to be stable. Most runs were approaching 

Sead, Stauc conditions at 600 time steps except the flat 
bottom cases which were running out to almost 1800 time 
Steps before settling down into steady state conditions. 

It should be noted that running the model as a flat bottom 
basin under these conditions was simulating the dynamics of 
Munk's model (1963). Therefore the results were somewhat 
DURCH BLeC seme CneSe Yuns. 

vie next set of runs were made with only source-sink 


eee necot Merons ON the western boundary. These runs 
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were intended to determine two items. The first was a 
Scivortcoonect the pattern of the resultant flow to that 
caused by a constant torque driving force. These results 
were compared for the flat bottom and the 20% bathymetry 
cases. Secondly, these runs were designed to test the 
model once again for stability but using a different 
Pouclie suet Lon ane the tlat bottom, 208 and 50% bathymetry 
cases. 

Waeer these eonditions, there appeared to be much less 
bathymetric dependence among the different runs. The first 
run was with flat-bottom conditions. As would be expected, 
With the S0s@bathymetry example, the transient topographic 
Rossby waves (Veronis 1966) had relatively high velocities 
ame the carculation pattern settled down wery austere in 
approximately 200 time steps. However, all cases had 
reached a steady-state condition by 600 time steps. 
beth the present conditions, the model appeared to be 
He telenulyeesensitive to numerical and configurational 
constraints while remaining stable throughout the timed 
runs. 

Before going on to the runs using the climatological 
data with the addition of the other variables, it was felt 
necessary to make further checks on the circulation patterns 
developed. “This was done to check consistency with the 
expected results due to the known dynamics involved to 
eisepewles “Ihe first case rum was Conducted with a 


Sonsteanetouque apolicd as the driving function with 


De 





no sources and sinks added. Non-linearities and bottom 
itieetTonmivercelett out (Gilgs="@5 = 0). Thé®*run started 
with the water initially at rest. A constant depth con- 
dition was imposed while the coriolis parameter varied 
Teeeldureuide. Under these conditions, a strong gyre 
within the basin was expected with some intensification 
five WeStern portion similar to Munk s moce 1.11963 )% 
Figure 6 shows that this was basically what was achieved. 
The steady state condition of this flow represents a balance 
Between the torque: applied by the wind and the dissipation 
of same by the lateral friction working on the model, 

The second case to consider is the same as the first 
fei wnemexception that 20, of the bathymetric deviations 
from the mean depth were allowed to interact. In this case, 
it was expected that the flow pattern would still basically 
be cyclonic although now, topographic Rossby waves were 
possible (Veronis 1966). Under these conditions, the flow 
could be expected to be broken into sub gyres centered in 
the basins around the Shirshov and Bowers Ridges. The 
results showed all this (Poem enmbirine tine Spin-up, 
propogation of a Rossby wave was noticed moving across the 
Study area from east to west. Cross ridge flow was observed 
in this "eases 

The third case involved increasing the bathymetry to 
Uo cCHpeetcCmmmmnnS did mot change the pattern drastically. 
However, the spin-up time to reach steady state conditions 
was reduced as could be anticipated and less transport was 


MoOLivced . 
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MiemitestmscerlesS USCd SOUrces and sinks on the western 
boundary as forcing functions. The wind torque was taken as 
zero. Again, non-linearities and bottom friction were ex- 
SdCimmincmilrSt Cases looked at a flat bottom situation. 
iiem low vwasedraven Dy tite input-output, source-sink initial 
conditions. Two small gyre patterns were established along 
the western boundary of the area (Figure 8). The gyre 
patterns along the boundary were a result of the localized 
PoUmec-oljmednivane toree. The gyre in the northwest corner 
WocederecouUlenor the erfects of the driving function combined 
With the geometry of the model. 

iiWjemseceondeease Of this Scries looked at the basic 
Sweuartoneswath 20% of the bathymetry included (Figure 9). 
The gyres created were the same as in the flat bottom case, 
with some noticeable effect from the bathymetry. The obvious 
SencctS un the bathymetry involve the magnitude of the 
Praisponresm AS expected, the transport volumes within the 
Majer portwen of the basin were reduced. The large-scale 
aneleyCcloniG eyre around hicemiagonmeporcion Of the area was 
a result of the bathymetry forcing the basic flow across the 
ines along the boundary. The small gyre in the south is 
being contained there by the Bowers Ridge. 

In comparing the effect of wind and of sources and sinks 
separately at this point, it appeared that the major driving 
ferece for the medel was the wind as opposed to the sources 


and sinks which appeared to have large effects only within 


Gilsmeelocalized a1rcawmeeroth of these series were now rerun 
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extending each case to 1800 time steps (900 days) in order 
POmoDoeivemrne Long-term elfects of the two driving forces. 
ieecatemcases, Circulation patterns remained qualitatively 
iNcoascmMemmO@lece tne Steady state conditions were met. The 
Berlowine section wili discuss the results of the next series 
©f runs. These were conducted with the climatological wind 
aeaeeas Grscussca Carlier, applying the stress to the 
surface along with the source-sink driving function. 
Comparison runs were made with and without non-linearities 


al sbDOwtom friction. 





IV. ACTUAL WIND DATA RESULTS 


iieeomsetics COlatsted of nine runs. They were con- 
ducted to investigate the contrast between a July and a 
Janwanry wind field driving force and the effects of varying 
bottom friction and non-linearities. Sources and sinks 
were added during the series for further comparisons. 

Wiest @mcote Lest Was Conducted with two runs. Both 
runs used July climatological wind data as the sole driving 
function. Non-linearities and bottom friction were allowed 
Dommecticine2enomme ine COrlolis parameter in ali cases varied 
With latitude as in all previous cases. The first run was 
a flat bottom case and the second run was with 20% bathymetry. 
In the first run, there was an intensified flow noticeable 
in the northwest corner of the area, as in the exploratory 
runs. This indicated that wind was forcing a westward 
intensification which was to be expected with the predominant 
Peevallingesouchwesterly winds Over the basin. The wind 
PrclicMmnacdme aim anrerey eG LOnle perturpation in the southwestern 
corner. This introduced some perturbations in the flow 
iveld anemevemmpboduced a negative gyre in the southwestern 
corner (Figure 10). After 600 time-steps, this case had 
reached an approximately steady state. 


The second run reached steady state conditions around 


420 time-steps (210 days). With the bathymetric features 


applied in this case, it was evident that the flow was 
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Dawe meury-cependent (Fipure 11). The gyre noticed in the 
first case moved around somewhat to align itself with the 
depth contours. Another anticyclonic gyre appeared to the 
east on the other side of the Bowers Ridge. The general flow 
PieemimioncyClOnic With perturbations throughout, evidently 
POMmclemeawenymMetric Crtects,. Transport values are reduced 
meaone-tamten trom the first case. There was some indication 
Sreeross ridge flow over the Shirshov and Bowers Ridges. 
iiiemle«e set Comsisted of three runs, These all 
PoecedsoOnty July wind field data and excluded non-linearities. 
itemourpescwor these runs Was to experiment with different 
mune ss ot bottom £riction. Including bottom friction in the 


equations trequires an equation of the form: 


of -_y¢ (IH 


where: 


This, in turn, requires an exponential solution of the form: 


K 4 


€=E 5° (1 2] 
which is a time dependent function. This coefficient can 
then physically be related to a l/e decay in the vorticity 
piela, Slimerctore values were determined for 1 (ene K a 
where K = C3 - Rand h was taken as the mulansingren dient hn sot atic 


model, 30 meters. Three cases were investigated where t was 


Peelemento so, 10, and 20 days. These corresponded to @ 
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decay time Within the model of 10, 20, and 40 time steps 
mespeectrye ly. 

All three cases reached steady state conditions between 
Wie idnosomemre Stevsn(150-175 days). In all cases the 
flow patterns were very similar indicating, as expected, that 
the bottom friction, once applied does not quantitatively 
eel tiemilow aS 1t 1S varied. In this set the flow did 
tend to follow the topographic aeidee aoe Giosely han 
Seserved previously. The value of C35 corresponding to a 
10 day period was finally chosen as the representative 
bottom LYriction parameter because the implied time scale 
seemed reasonable under the imposed conditions (Figure 12). 

livemmex: Set consisted of two runs. This set was 
Someneccamuec Imvesticate the effects of sources and sinks. 
The basic data used on each of the runs was the July wind 
mela, Ui wbathnymetry, CS for 10 days as determined in the 
Pots ScCtemeorlOlisS parameter varying with latitude and 
Somrec-sSummamadud £Or the western boundary. Initial conditions 
Pomel Sescumlesewere taken £rOm alprevious Series run with 
Gey es COULeeesimkwedata tfor 900 days in order to exclude the 
source-sink driven transients as aucomemrGating factor. 
(ne hiustmeumn Was Completed excluding non-linearities 
(Figure 13). In the second run, non-linearities were 
included |C |= —--*—|(Figure 14). 

ASDF 

Mic couemmiowepatterns showed gyres in the local- 

iecimeneamunere themsource-Sink cross-boundary transports 


pewemopelteGe nas Was tombe expected. An intense pattern was 
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set up at the interface between the motions due principally 
Pomemcmoolrces and sinks and those due to the winds. In 
POWeasese a esigniticane Lact is that it appeared that the 
SOurces and sanks were the dominant forces in the flow field. 
Pi tenwecOleideracion Of the relative driving forces, this was 
not totally unexpected. A negative gyre appeared northwest 
of Shirshov Ridge and the gyres previously observed on either 
side of the Bowers Ridge remained. Cross ridge flow was 
SvidemesimebOtmecases. In the second case, as was anticipated, 
the addition of the non-linearities smoothed the flow somewhat. 
litemiasteset also consisted of two runs. This time the 
January wind field was used. The first run was without 
non-linearities (Figure 15), and the second run included 
them (Figure 16). Smoothing was evident in the second run. 
Wanwaryetas the strongest wind fieldof the existing data and 
as such it developed stronger gyres throughout the basin 
area than the July wind. In this set, it became evident that 
the stronger wind field was influencing somewhat the force 
of the sources and sinks. The gyres along the western bound- 
apy were still present, as was expected in the local area of 
Piers OUnGecw ana soils oUt thie Strong@eanticyclonic gyre 
over the remainder of the Bering Sea basin had been reduced 
cide NOneeie weds ECa ne pOmnlonemim particular, had even been 
Pm placcduy ad cyelonieupertirbation. This perturbation was 
viewed as a direct result of the stronger wind field. In 
ae wo WmENeenOnoencarities appeared to have reduced the 


magnitude of the gyres along the Aleutian chain and had 
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Vemiiicely smoothed out the perturbation field as a result 
of the interaction of the true driving forces. This became 
noticeable along the eastern side of the basin where there 
Mewrex1Stea a possibility of a counter current due to the 
smoothing effect. 

The counter current hinted at has also been discussed 
by Coachman (1973) as a possibility during certain times of 
Mie year. ©it 15 not clear that the counter current in the 
model is dynamically similar to the one he described. 
However, in order to look at this possibility more closely, 
the next section investigated the results of adding free 


SMpmconaaterois along the eastern boundary. 


4§ 





Nome NvEoT LIGATION OF FREE SLIP BOUNDARY CONDITIONS 


The model had been run with only no-slip conditions 
imposed on all boundaries at this stage. This next series 
Mosmeenauceed Co investigate the results of adding free 
slip COnaERE ors along the open boundary on the eastern 
edge of the model. This area, in the the Bering Sea, is 
ali open water along thescontinental shelf and, for the 
Mest part, the depth of water is less than 300 meters. 

Two separate runs were made, one with and one 
Without bottom friction. The January wind field data 
was used along with the sources and sinks as the driving 
iOomcewmene Lun Without bottom Eriction (C3 = 0), Figure 17, 
Sewed smemose resemblance to Pighre 16 everywhere except 
athe ated Or direct interest. Alomg the eastern and 
northern edges of the model, a definite cyclonic gyre was 
forming. It was confined to the continental shelf region 
and the Shirshov Ridge area Ones ncoweoOuntLer Gurrent 
possibility still existed along the continental slope and 
shelf areas, primarily in water less than 1000 meters deep. 

In the second run, when bottom friction was added in, 
ine CyCtlenic yre Spreadsout more in the northern sector 
and was forced into the western area by the Shirshov and 
Bowers Ridges. Along the eastern edge, there appeared to 


be an attempt for the counter current to follow the 
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woorrepntec relict along the 1000 meter contour. This 
Muteecdrit right along the continental slope (Figure 18). 
Deen these runs ha@ non-litearitires ing@luded in the 
problem. The addition of free slip conditions along the 
eastern open boundary did not seem to make any appreciable 
difference in the flow pattern or in the magnitude of the 


transports. 
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Yi. CONGEUSTONS 


Following the setting up of the model, testing for 
sensitivity and stability were conducted to ensure the model 
Tevet kinegsproperiy., lo dem@this, an initial run under 
flat bottom conditions, with a constant wind force applied, 
was made. The results were not unexpected. The flow was 
cyclonic and had an obvious intensification along the 
western boundary similar to Munk's closed basin studies. 

Actual wind data was used to drive the model in the 
last series of runs. This was data compiled by month. 

July data was the weakest field and January was the strong- 
est. these were combined with source-sink transport exchange 
values along the western boundary. Comparisons were made 
between results of the July and January wind fields, 

while including and excluding bottom friction and non-linearities. 

In all cases when the bathymetric factors were added, 
the lew wceeme very dependent om the topographic relief. 

The Shirshov and Bowers Ridges tend to divide the area into 
three basins, the Kamchatka Basin, the Aleutian Basin, and 
an area between the two ridges. The overall circulation 
Pactern was civided into gyres, one in each basin. 

Initially, when looking at the wind and sources and sinks 
as separate driving forces, it appeared that the wind was 
the primary driving force. This also agreed with the results 


from previous researchers. The results of the Sorce-sink 
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driving forces alone was two principal cyclonic gyres 

along the western boundary and one large anticyclonic 

eee cnever the mwemainder of the basin. The wind force 
Paocuecd aclargcewcyclonic gyre over the entire basin. 

Once the two driving forces were put together, the apparent 
dominant force became the sources and sinks. 

In the case where the wind and Somece and sinks were 
used to drive the model, the flow became generally anti- 
cyclonic. There were two rather separate patterns, however. 
Along the western boundary the flow was definitely driven 
Date source-Sink input. As the flow got further from the 
western boundary, there was an area of intensification which 
Was a result of the interaction of the wind influence and 
piete ot the wseurces and sinks. Beyond that, came the major 
Ceaec tidtecomered the remainder of the deep sea basin area. 
This was an anticyclonic gyre with some perturbations in 
PicmmlClidweminese were viewed as@peanpe influenced” by the 
wind field. 

July wind data was used at first and this resulted in 
the pattern just described. These runs had bottom friction 
ain the €quations. The effects caused by the bathymetry were 
much larger than the effects observed by the bottom friction. 
Some of the runs had non-linearities included. These showed 
Piaiichiiiyea smOOthingsot the flow pattern. When the 
January Wind field data was used, the flow in the eastern 
portion of the area started to show the effects of the 
Saronc wind tield overriding the influence of the sources and 


Simks. 
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Tie timal set of runs were conducted to explore the 
differences between free slip and no slip conditions along 
the eastern boundary. With the stronger January wind 
applied in the model, the cyclonic gyre that appeared in the 
eastern portion also introduced the possibility of a 
counter current along the continental slope. By changing 
the initial conditions to allow free slip Hie thateawea a tne 
Goumner CcUnment tended to follow bathymetry a little closer 
but no more than when non-linearities were added. 

Highes spremescd Circlubatvom seneme (figure 19) was 
feereloved for the summer months and as such, the July 
Soneostue Huns were used for comparison. Generally, there 
Mmemerareas Or agreement as well as disapreement. Along 
the western boundary, the flow pattern agree quite well. 
In the northwestern corner of the Kamchatka Basin, both 
Patterns Showed an anticyclonic gyre. the two patterns 
disagreed in the overall flow. Hughes shows a cyclonic 
tendency where the model has an anticyclonic circulation; 
Pemercutaning, §j;USt west of the Shirshov Ridge where the 
model showed a large negative gyre approaching two Sverdrups 
me the summer and three to four Sverdrups in the winter. 
Evident in both instances, however, is the fact that gyres 
throughout the Bering Sea deep water basin are predominantly 
dependent upon topographic relief. Cross ridge flow is 
also evident in both. 

ihesie comparisons have shown that for a thorough 


understanding of the dynamics of the Bering Sea, study must 
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Gonmuinue along both these lines. More direct measurements 
mbmewieceessary, particularly along the castern boundary at 
the edge of the continental Slope These would aid in prob- 
ing the counter current suspected in that area. In the 
numerical modcling, studies are presently continuing on 

a time-dependent wind driven model. This Study uses twelve 
months of climatological data. Further studies should be 
conducted with investigation into source-sink data along 


maer other boundaries of the model. 
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